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ABSTRACT 

Context. The magnetic field permeating the solar atmosphere is generally thought to provide the energy for much of the activity 
seen in the solar corona, such as flares, coronal mass ejections (CMEs), etc. To overcome the unavailability of coronal magnetic 
field measurements, photospheric magnetic field vector data can be used to reconstruct the coronal field. Currently, there are several 
modelling techniques being used to calculate three-dimensional field lines into the solar atmosphere. 

Aims. For the first time, synoptic maps of a photospheric-vector magnetic field synthesized from the vector spectromagnetograph 
(VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) are used to model the coronal magnetic field and estimate 
free magnetic energy in the global scale. The free energy (i.e., the energy in excess of the potential field energy) is one of the main 
indicators used in space weather forecasts to predict the eruptivity of active regions. 

Methods. We solve the nonlinear force-free field equations using an optimization principle in spherical geometry. The resulting three- 
dimensional magnetic fields are used to estimate the magnetic free energy content E free = E nm - E pot , which is the difference of the 
magnetic energies between the nonpotential field and the potential field in the global solar corona. For comparison, we overlay the 
extrapolated magnetic field lines with the extreme ultraviolet (EUV) observations by the atmospheric imaging assembly (AIA) on 
board the Solar Dynamics Observatory (SDO). 

Results. For a single Carrington rotation 2121, we find that the global nonlinear force-free field (NLFFF) magnetic energy density 
is 10.3% higher than the potential one. Most of this free energy is located in active regions. 

Key words, magnetic fields - Sun: corona - Sun: photosphere 


1. Introduction 

In the solar coronal plasma, magnetic energy is the prime energy 
reservoir that fuels the dynamical evolution of eruptive events, 
but it remains an open question to how the magnetic energy is 
released. The amount of energy associated with the magnetic 
field is much larger than other energy sources, and the dynamics 
of the coronal configuration is determined by the evolution of its 
magnetic field (Forbes 2000; Low 2001). 

Free energy is defined as excess of energy as compared 
with potential field. One can, in principle, release more energy 
than free energy by, for example, annihilating magnetic field. 
Numerically, one can also get pre-post flare difference in energy 
larger than free magnetic energy if the magnetic field at the pho- 
tosphere changes (e.g., some flux elements disappear). It quan- 
tifies the energy deviation of the coronal magnetic field from its 
potential state (Metcalf et al. 2005; Regnier 2007; Aschwanden 
2012). The magnetic free energy is stored in the form of electric 
currents flowing along the magnetic field. Free magnetic energy 
of solar magnetic fields can be affected by several processes such 
as photospheric shearing flows, magnetic flux emergence and 
magnetic reconnection (e.g., Welsch 2006; Fang et al. 2012). 

To understand the role that the magnetic field plays in ener- 
gizing the solar corona, it is important to calculate the amount 


of free energy needed to quantify the energy budget in a catas- 
trophic energy release event, as well as for estimating upper lim- 
its in forecasting individual events in real-time. Using various 
extrapolation techniques for the coronal magnetic field under the 
assumption of force-free fields, the spatial and temporal evolu- 
tion of the coronal magnetic free energy during solar flares has 
been extensively studied (Regnier & Canfield 2006; Guo et al. 
2008; Jing et al. 2010; Tadesse et al. 2012a; Meyer et al. 2013). 
From the 3D coronal magnetic configurations, we can derive the 
magnetic energy in the corona as 

E m = J- j B B r 2 sin#drd#d0. (1) 

Jv 

The free magnetic energy in spherical geometry is calculated by 
computing the nonpotential field B n mf(r, 6, 0) with a numerical 
nonlinear force-free field (NLFFF) code and a potential field 
Bp 0i (r, 6, 0) for the same photospheric boundary data, so that the 
difference of the magnetic field energy density in the coronal 
volume V encompassing the active regions of interest can be 
quantified as E Ufm = E Unm - E Mpot . 

Since the corona is optically thin, direct measurements of 
the 3D magnetic field are very difficult to implement and inter- 
pret. Therefore, the present observations for the magnetic fields 
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based on the spectropolarimetric method (the Zeeman and the 
Hanle effects) are limited to low layers of solar atmosphere (pho- 
tosphere and chromosphere). Even if direct measurement tech- 
niques for the 3D magnetic fields in the chromosphere and the 
corona have considerably improved in recent decades (Lin et al. 
2000, 2004; Liu & Lin 2000), further developments are needed 
before accurate data are routinely available. Thus, the problem of 
measuring the coronal field and its embedded electrical currents 
leads us to use numerical modelling to infer the field strength 
in the higher layers of the solar atmosphere from the measured 
photospheric field. 

Lorce-free extrapolation of photospheric magnetic fields is 
currently used as the primary tool for the modelling of coro- 
nal magnetic fields. In this model assumption, the corona mag- 
netic forces are dominant so that all nonmagnetic forces, such 
as the plasma pressure gradient and gravity, can be neglected in 
the lowest order (Gary 2001). This implies that, if appreciable 
currents are present, these must be aligned with the magnetic 
field, since the resulting Lorentz forces could otherwise not be 
balanced by the nonmagnetic forces. The equilibrium structure 
of the static coronal magnetic field can be described using the 
force-free assumption as 

(V xB) = aB ^ (V xB)xB = 0 (2) 

V • B = 0, (3) 

where B is the magnetic field. The force-free parameter a of 
Eq. (2) can be a function of position, but the combination of 
Eqs. (2) and (3) ( B • Vcr = 0) requires that a be constant along 
a given field line. Potential ( a = 0) and linear force-free fields 
(whenever a is constant everywhere in the volume under con- 
sideration) can be used as a first step to model the general struc- 
ture of magnetic fields in the solar corona. Practically, the pre- 
emptive magnetic fields are nonlinear force-free fields ( a being 
a function of position) as supported by both observational and 
theoretical reasons. Lor details of those models, we direct the 
readers to Wiegelmann & Sakurai (2012). 

Nonlinear force-free field codes have been routinely ap- 
plied to the reconstruction of the coronal field of a single active 
region using the Cartesian geometry. In that case, the curva- 
ture of the solar surface does not play an important role. The 
Solar Dynamics Observatory (SDO) mission has made repeated 
observations of large scale events in which connections be- 
tween widely separated active regions may play fundamental 
role (Martens et al. 2012). Even before SDO, it was known that 
large-scale connectivity is important for solar eruptive and non- 
eruptive activity (e.g., studies of sympathetic flares, transequa- 
torial loops (Pevtsov 2000), effects of distant active regions on 
large-scale coronal brightness (Pevtsov & Acton 2001), eruption 
of filaments triggered by remote flux emergence, and evolution 
(e.g., Balasubramaniam et al. 201 1). Therefore, these needs mo- 
tivate us to implement a NLFFF procedure in spherical geometry 
(Wiegelmann 2007; Tadesse et al. 2009, 2011, 2012a, b, 2013b, 
2014; Guo et al. 2012; Amari et al. 2013) 

Nonlinear force-free reconstruction of the global solar mag- 
netic field from line-of- sight component of synoptic magne- 
tograms has been developed by Contopoulos et al. (2011). In this 
study, we estimate the free magnetic energy for global corona 
using the synoptic vector magnetograms from SOLIS/VSM. We 
compare the extrapolated potential and NLFFF magnetic loops 
with extreme ultraviolet (EUV) observations by the atmospheric 
imaging assembly (AIA) on board SDO. This comparison helps 
to identify whether the NLFFF model reconstructs the magnetic 
configuration better than the potential field model in the global 


scale. In this paper, we present some descriptions of the dataset 
used for analysis in Sect. 2. The spherical optimization proce- 
dure used for modeling 3D magnetic field in global corona is 
presented in Sect. 3. Then, we present results of our studies 
in Sect. 4. A summary and discussions are finally presented in 
Sect. 5. 


2. Instrumentation and data set 

In this paper, we use the first ever synoptic vector magne- 
tograms, which include Carrington maps of the three compo- 
nents of the magnetic field vector, the radial B r , the poloidal B e , 
and the toroidal B to first perform a global nonlinear 
force free field extrapolation based on optimization algorithm 
(Wiegelmann 2007; Tadesse et al. 2009). The synoptic vector 
field maps are derived using daily full disk photospheric vec- 
tor magnetograms from vector spectromagnetograph (VSM) in- 
strument on the synoptic optical long-term investigations of the 
Sun (SOLIS), a synoptic observing facility (Keller et al. 2003; 
Balasubramaniam & Pevtsov 2011). More detail about these 
synoptic maps and their properties can be found in Gosain et al. 
(2013). Here, we briefly describe the instrument characteristics 
and the full disk vector field observations. 

The vector spectromagnetograph routinely obtains full disk 
magnetograms in photospheric and chromospheric lines as a 
part of the synoptic program of National Solar Observatory, 
known as the NS O Integrated Synoptic Program (NISP). To ob- 
tain full disk photospheric vector magnetograms, SOLIS/VSM 
measures Stokes ( S = /, Q, U , V) profiles in photospheric 
Fe I 630.15-630.25 nm line pair. A single full disk scan 
(2048 scan lines) takes only about 20 minutes, thanks to the 
long slit of the spectrograph, which intersects the solar disk 
from one limb to another in one shot. The spatial sampling 
is 1 arcsec per pixel with square pixels. The spectral sampling 
is 2.4 pm per pixel. A single Stokes set (/, Q, U, V ) per slit po- 
sition is obtained typically in about 0.6 s. The telescope itself 
is designed to be free of instrumental polarization by employ- 
ing symmetric optical configuration and performing polarization 
modulation just after the prime focus after the slit. Dual beam 
analysis of polarization using a polarizing beam splitter avoids 
seeing that is induced cross talks in the signal. The polarime- 
ter calibration is done routinely to calibrate the Stokes vector 
for cross-talks. The signal- to-noise ratio in the continuum of 
the Stokes profiles is typically >1000. The magnetic field vec- 
tor is inferred from the calibrated Stokes profiles by performing 
inversion in the framework of the Milne-Eddington model for 
stellar atmosphere, which follows Unno-Rachkovsky formalism 
(Skumanich & Lites 1987). 

Only pixels with a polarization signal above the thresh- 
old of 0.1% of continuum intensity, Ic, are inverted to obtain 
the magnetic (field strength, inclination angle, and azimuth an- 
gle) and thermodynamic (e.g., Doppler width, Doppler velocity, 
source function, temperature) parameters. The threshold of 0. 1 % 
of Ic corresponds to typical noise level in the continuum. Using 
this threshold avoids fitting profiles buried in the noise. 

Further details about instrument and pipeline reduction steps 
can be found elsewhere (e.g., Jones et al. 2002; Henney et al. 
2006; Balasubramaniam & Pevtsov 2011). The noise in SOLIS 
magnetograms is estimated to be a few gauss in longitudinal 
and 70 G in the transverse field measurements (Tadesse et al. 
2013a). The 180 degree azimuth ambiguity is resolved using 
a different (faster) ambiguity resolution method developed by 
Rudenko & Anfinogentov (2011). 
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3. Magnetic field modelling 

For modelling the coronal magnetic field in a global scale, we 
use the variation principle originally proposed by Wheatland 
et al. (2000). Wiegelmann (2007) has developed an optimiza- 
tion method to reconstruct the NLFFF for global solar corona 
by minimizing an objective functional L that combines Lorentz 
forces and the divergence of the magnetic field in spherical ge- 
ometry. The code has been tested with semi-analytic force-free 
solutions (Low & Lou 1990). If the functional is minimized to 
zero, Eqs. (2) and (3) are satisfied simultaneously. The optimiza- 
tion procedure in the spherical geometry has been implemented 
by Tadesse et al. (2009) for a restricted area with large field of 
views. Later, Tadesse et al. (2011) modified the objective func- 
tion of the optimization method (Wiegelmann & Inhester 2010) 
for spherical geometry as 

L — Lf + Ld + vLphoto (4) 

u = I s -2 |(V x B) x B| V sin Odrdedcp 
Jv 

Ld = I |V • B^r 2 sinOdrdOdfi 
Jv 

iphoto = J' {B - B obs ) • W(0, 0) • (B - B obs )r 2 sin 6d6d<p, 

where Lf and Ld measure how well the force-free Eqs. (2) and 
divergence-free (3) conditions are fulfilled, respectively. The 
main reason for modification of the code was that we need to 
deal with boundary data of different noise levels and qualities, or 
even miss some data points completely. Hence, the third integral, 
L photo , is the surface integral over the photosphere, which allows 
us to relax the field on the photosphere towards a force-free so- 
lution without too much deviation from the original surface field 
data. 

The SOLIS/VSM provides full disk vector magnetograms 
from which synoptic Carrington maps of the vector magnetic 
field components are synthesized. However, for pixels below the 
polarization threshold the inversion was not performed and field 
data will be missing for these pixels (see Fig. 1). Typically, the 
field is missing where its magnitude is small; thus, these pix- 
els would have a small impact on the model even if they were 
measured correctly. Within the error margin of a measured field 
value, any value is just as good as any other, and from this range 
of values, we take the value that fits the force-free field best. To 
treat those pixels with missing data, we used the diagonal matrix, 
W(0, 0), which gives different weights to the observed surface 
field components depending on the relative accuracy in measure- 
ment. We choose W = (B T / ma x(B T )) 2 . This seems to be a rea- 
sonable choice, as the measurement error in the transverse field 
is higher in weak field regions (Wiegelmann & Inhester 2010). 
In this sense, the missing data is considered the most inaccu- 
rate and is taken into account by setting W(0, 0) to zero in all 
elements of the matrix. 

The photospheric magnetic field has a plasma-/? of order 
unity, which does not satisfy the force-free condition (Gary 
2001). Therefore, the vector magnetogram data are inconsis- 
tent with the force-free assumption, which is absolutely essential 
condition for NLFFF extrapolation. To find suitable boundary 
conditions for the NLFFF field modelling, we have to prepro- 
cess the measured synoptic vector magnetograms by using a pre- 
processing scheme developed by Tadesse et al. (2009) in spher- 
ical geometry. This preprocessing scheme removes forces and 
torques from the boundary and approximates the photospheric 


magnetic field to the low plasma-/? force-free chromosphere. For 
a detailed description of the current code implementation, we 
refer to Wiegelmann (2007) and Tadesse et al. (2011). Tadesse 
et al. (2012a) has compared the difference between the boundary 
data before and after the preprocessing and found that the maxi- 
mum change in the transverse field is 30 G for SOLIS/VSM data. 
The final preprocessed field values on the boundary are consis- 
tent with the observed data within noise levels (noise due to fluc- 
tuations in intensity), which are about 1 G and 70 G for longitu- 
dinal and transverse components, respectively. 

4. Results 

This study requires extrapolating the three-dimensional poten- 
tial and NLFFF coronal fields from the photospheric boundary 
on global scale. We use synoptic maps of photospheric vec- 
tor magnetic field observed during 4-31 March 2012. During 
this observation, there were about 25 active regions all over the 
solar globe. To use our spherical optimization code for global 
corona, we adopt a uniform spherical grid using r, 6, and 0 with 
n r = 300, ng = 450, and n $ = 900 grid points in the direction 
of radius, latitude, and longitude, respectively, with the field of 
view Of \r m \ n = 1 Rq r max = 2.5 Rq] X [0 m in = — 90 ! ^max = 
90°] x [0 min = 0° : 0 max = 360°]. To avoid the mathematical 
singularities at the poles, we do not use grid points exactly at 
the south and north pole but set them half a grid point apart at 
#min = -90° + d6 /I and 6 max = 90° - d6/2 (see, Wiegelmann 
2007). We set W(0, 0) = 0 for those pixels where there are no 
measurements on the polar regions. The code solves the NLFFF 
equations in the bounded domain between 1 R Q and the source 
surface at 2.5 R Q . The outer boundary is kept fixed using the ini- 
tial potential field values. All current-carrying field lines have 
to be confined the volume. The domain outside 2.5 R Q is not 
included in the model, because the force-free approach is no 
longer justified here. The magnetic field extrapolation has been 
carried out by excluding the polar regions for almost global 
corona where the magnetic field measurements are highly being 
influenced by noise. Before performing NLFFF extrapolations, 
we use the preprocessed radial magnetic field component B r 
to compute the corresponding potential field using a spherical 
harmonic expansion for initializing our spherical NLFFF code 
during relaxation towards a force-free state in the computational 
volume. 

The main purpose of this work is to study the structures 
of the global potential and NLFFF magnetic fields and to es- 
timate free magnetic energy available to power solar eruptions 
during Carrington rotation 2121. In addition, we compare which 
of those two models best agrees with observation in global en- 
vironment. To do this, we plot the selected fieldlines of the 
potential and NLFFF models in Fig. 2. We overlay the field 
lines with an AIA 171 A and a composite (171 A, 193 A, 
and 211 A) image. The field lines of the potential and NLFFF 
models are reconstructed from the same footpoints. Some of 
potential field lines in Fig. 3a have an obvious deviation from 
the observed EUV loops. There are more field lines open to 
the upper boundary in NLFFF than potential fields. Most of 
AIA coronal loops are best overlaid by NLFFF lines than po- 
tential ones (see Fig. 3b). Therefore, the qualitative compari- 
son between the model magnetic field lines and the observed 
EUV loops indicates that the NLFFF model provides a more 
consistent field for global corona magnetic field reconstruction. 
Figure 4 shows that there is over all similarity between the 
corresponding NLFFF model field lines and image of the sun 
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B(0) CR-2121 



0 100 200 300 


Fig. 1 . Synoptic Carrington maps of the vector magnetic field components are shown for CR-2121. The panels from top to bottom show the 
distribution of the B r , B $ and B e components, respectively. The B r map is scaled between ±100 G, and the B $ and B e maps are scaled to +20 G. 
The positive values of B r , B $ and B e point, respectively, upward, westward and southward. 
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(a) AIA 171A 



(b) AIA 171, 193, and 211 A composite 



(c) Potential field over AIA 171 A 



(d) Potential field over AIA composite 



(e) NLFF field over AIA 171A 



(f) NLFF field over AIA composite 

Fig. 2. Global field lines of the potential field model and the NLFFF model overlaid on the AIA 171 A and composite AIA (171, 193, and 21 1 A) 
images. Green and red lines represent open and closed magnetic field lines, respectively. 


observed by the SOHO/LAS CO C2 coronagraph. However, one 
can see that NLFFF does not well represent a linear structure 
above a coronal helmet in the low-left corner in Fig. 4b. This 
could be due to the missing data from the polar region. 


In addition to the above, to compare the degree of disagree- 
ment between the two model vector field solutions in the global 
corona volume that are specified on the identical sets of grid 
points, we use the vector correlation metric (C vec X which is also 
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Fig. 3. Field lines of a) the potential field model and b) the NLFFF model around ARs 1 1429 and 1 1 430 overlaid on the AIA 171 A image. Green 
and red lines represent open and closed magnetic field lines, respectively. 



(a) (b) 


Fig. 4. Magnetic field line skeletons a) of the entire solar globe from the NLFFF model and b) image of the sun observed by SOHO/LASCO 
C2 coronagraph at 16:33UT. 


used to the standard correlation coefficient for scalar functions. 
The correlation was calculated (Schrijver et al. 2006) from 


C V ec — 


E; Vi ■ Ui 


(ZikPZiN 2 ) 


1 / 2 ’ 


( 5 ) 


where Vi and m* are the vectors at each grid point i. If the vector 
fields are identical, then C vec = 1; if 0 / -L M; , then C vec = 0. The 
degree of convergence towards a force-free and divergence-free 
model solution can be quantified by the integral measures of the 
Lorentz force and the divergence terms in the minimization func- 
tional in Eq. (4), which are computed over the entire solar globe. 
The Lf and L & of Eq. (4) measure how well the force-free and 
divergence-free conditions are fulfilled, respectively. In Table 1, 
we provide some quantitative measures to rate the quality of 
our reconstruction. Column 1 names the corresponding models. 
Columns 2-3 show how well the force-balance and solenoidal 


Table 1. Evaluation of the reconstruction quality for the potential field 
and NLFFF models. 


Model 

U 

L d 

Tphoto 

c 

'-'vec 

Potential 

0.000 

0.000 

0.001 

1 

NLFFF 

0.391 

0.697 

0.302 

0.893 


Notes. We have used spherical grids of 300 x 450 x 900. 


conditions are fulfilled for both models. Figure 5 shows how 
well the functional L converge to zero during iteration process. 
In the last column, the vector correlation shows that there is dis- 
agreement between the two model field solutions. 

The energy stored in the magnetic field as a result of a field 
line stressing into a nonpotential configuration has been identi- 
fied as the source of flare energy. Therefore, to understand the 
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Fig. 5. Evolution of the entire functional L (blue line) and its three terms 
in Eq. (4) during the optimization process. The black line corresponds 
to Lp h oto, the red line to L f , and the green line to L d . 


Table 2. Magnetic energy associated with the 3D potential and NLFFF 
field configurations calculated from synoptic vector magnetogram. 


Model 

F to tai(10 33 erg) 

Ff r ee(10 33 erg) 

Potential 

41.6 

0 

NLFFF 

46.4 

4.8 


physics of solar flares, including the local reorganisation of the 
magnetic field and the acceleration of energetic particles, we 
have to first estimate the free magnetic energy available for such 
phenomena. This free magnetic energy can be converted into ki- 
netic and thermal energy. We estimate the free magnetic energy, 
which is the difference between the extrapolated NLFFF and the 
potential field with the same normal boundary conditions in the 
photosphere. We therefore estimate the upper limit to the free 
magnetic energy associated with the coronal currents of the form 

Efree = ££ - S 2 pot )dV. (6) 

The parameters 5p 0 t and represent the potential and 

NLFFF magnetic field, respectively. Our result for the estimation 
of free-magnetic energy in Table 2 shows that the NLFFF model 
has 10.3% more energy than the corresponding potential field 
model. Figure 6a shows a synoptic map of the global free mag- 
netic energy relative to the potential energy on the photosphere. 
In Fig. 6b, we plot the free magnetic energy density on the lat- 
itudinal plane 6 = 20° relative to potential field in the coronal 
height. There are strong free energy concentrations above each 
active region over the solar globe. 

The use of preprocessed boundary data and measurement 
error due to noise may influence the estimation of free mag- 
netic energy (Tadesse et al. 2012a). As the preprocessing proce- 
dure filters out, small-scale surface field fluctuations, the mag- 
netic energy of the NLFFF obtained from the preprocessed 
boundary data is smaller than the corresponding energy with- 
out preprocessing. The magnetic energy computed from the 
original data without preprocessing is 5 x 10 33 erg, which is 
about 4.2% higher than that obtained from the preprocessed 
boundary data. However, this energy does not correspond to the 
NLFF magnetic field solution, since the original boundary data 
without preprocessing is not a consistent boundary condition for 
NLFFF modelling. 



(a) Synoptic map of free magnetic energy 



0 100 200 300 

I_Ongitude[ Degree] 

(b) Free magnetic energy in height 


Fig. 6. a) Synoptic map of the global free magnetic energy relative to 
the potential energy and b) free magnetic energy map on the latitudinal 
plane 6 = 20° relative to potential field (increasing density from black 
to white). 


5. Conclusion and outlook 

Most of the NLFFF procedures are implemented in a Cartesian 
coordinates. Therefore, both potential and nonlinear force-free 
field (NLFFF) codes in Cartesian geometry are not well suited 
for larger domains, since the spherical nature of the solar surface 
cannot be neglected when the field of view is large. Therefore, 
it is necessary to implement a NLFFF procedure in spherical 
geometry for use when large-scale boundary data are in use. 

In this study, we have investigated the coronal magnetic 
field and free magnetic energy associated with the global corona 
by analyzing Carrington synoptic maps of the photospheric 
vector magnetic field, which are synthesized from the vector 
spectromagnetograph (VSM) on Synoptic Optical Long-term 
Investigations of the Sun (SOLIS). This has been observed dur- 
ing 4-31 March, 2012. The Carrington rotation number for this 
observation is 2121. During this particular observation, there 
were about ten active regions distributed across the globe. We 
have used our spherical NLFFF and potential codes to compute 
the magnetic field solutions over global corona. This is the first 
NLFFF magnetic field extrapolation for the global corona using 
real data. 

We have compared the magnetic field solutions from 
both potential and NLFFF models. The qualitative compari- 
son between the model magnetic field lines and the observed 
EUV loops indicates that the NLFFF model provides a more 
consistent field for global corona magnetic field reconstruction. 
For this particular Carrington rotation, we find that the global 
NLFFF magnetic energy density is 10.3% higher than the poten- 
tial one. For the future, we plan to study the evolution of global 
free magnetic energy. In this study, most of this free energy is 
located in active regions. 
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